Collective excitations and the nature of Mott transition in undoped gapped graphene 
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Particle-hole continuum (PHC) for massive Dirac fermions in presence of short range interactions, provides 
an unprecedented opportunity for formation of two collective split-off states, one in the singlet and the other 
in the triplet (spin-1) channel in undoped system. Both poles are close in energy and are separated from the 
continuum of free particle-hole excitations by an energy scale of the order of gap parameter A. They both 
disperse linearly with two different velocities reminiscent of spin-charge separation in Luttinger liquids. When 
the strength of Hubbard interactions is stronger than a critical value, the velocity of singlet excitation which 
we interpret as a charge boson composite becomes zero, and renders the system a Mott insulator. Beyond this 
critical point, the low-energy sector is left with a linearly dispersing triplet mode - a characteristic of a Mott 
insulator. The velocity of triplet mode at the Mott criticality is twice the velocity of underlying Dirac fermions. 
The phase transition line in the space of U and A is in qualitative agreement with a more involved dynamical 
mean field theory (DMFT) calculation. 

PACS numbers: 73.22.Pr, 73.20.Mf, 75.10.Kt 



INTRODUCTION 



Graphene is a huge network of sp^ bonds formed between 
carbon atoms on a two-dimensional (2D) sheet f\\\. The low- 
energy effective model for electrons in this system is the 
2h-1 dimensional Dirac theory which contains an energy scale 
proportional to vp - the Fermi velocity of electrons pi] . Ef- 
fects of interactions and/or disorder can be taken into account 
on top of this non-interacting fixed point Breaking the 
sub-lattice symmetry gives rise to a gap in the single-particle 
excitations spectrum |3|], which introduces the gap parame- 
ter A as another energy scale. The gapped graphene can be 
modeled with 2h-1 dimensional massive Dirac fermions fl- 
Gapped graphene supports interesting new class of excita- 
tions in the spectrum such as, domain walls IH. Moreover 
the single-particle gap (mass) gives rise to suppression of 
Coulomb fields at nano-meter length scales 161]. The gap 
could also give rise to universal nonlinear optical response ■ 
When the many-body interactions of various form are added 
to the massive Dirac theory, the nature of many-body exci- 
tations becomes even more interesting. Recent ab-initio es- 
timates of the strength of the Coulomb repulsion in various 
forms of graphene suggests that the Hubbard parameter in 
graphene can be quite remarkable which introduces a third 
energy scale, U. Therefore, an interesting question here would 
be the effects of local Hubbard type interactions on the excita- 
tion spectrum of gapped graphene, and the nature of transition 
to Mott insulating phase in gapped graphene. Our recent full- 
fledged DMFT investigation of the so called ionic-Hubbard 
model on the honeycomb lattice suggested three phases [jgt]: 
(i) Band insulating phase for A ^ U. (ii) Mott insulating 
phase for A <C ?7, and (iii) a semi-metallic phase for U ^ A. 
The above study is about the nature of ground state. In or- 
der to study the connection of excited states and the ground 
state, one notes that; deep in the Mott insulating phase, the 
energy scale U is expected to be dominant over A, such that 
the low-energy excitations are expected to be spin fluctuations 



arising from super-exchange interaction induced by the large 
U. We would like to focus on, the evolution of the collective 
spin and charge dynamics in this system as a function of the 
Hubbard U for a system with a non-zero gap parameter A. It 
turns out that the presence of a non-zero gap parameter facili- 
tates the separation of two collective states each of which has 
a different velocity, reminiscent of the spin-charge separation 
in Luttinger liquids. 



FORMULATION OF THE PROBLEM 

The Hamiltonian we consider here is the so called ionic- 
Hubbard model, which is defined by 



H = hvF ^ V'L ['^■^ + '^^^] V'ks' 
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where nj^ = fjsfjs "^^^^ I = (^jb corresponding to num- 
ber operator at the j'th unit cell at sub-lattices. A, B, respec- 
tively. The spinor notation tjjj^,, = (aj^^, 6j^^) has been used. 
Here s =t, i denotes the z component of the physical spin, 
k = {kx,ky) is a two-dimensional momentum vector, and 
a stands for Pauli matrices. U is the strength of the on-site 
Coulomb interaction known as Hubbard parameter, and A is 
the mass parameter. The quadratic part of this Hamiltonian 
can be diagonalized by a simple unitary transformation to the 
basis of conduction (+) and valence (— ) states, and the corre- 
sponding eigen-values are given by. 



(1) 



e± (k) = ±VF V|ftk|2 + A2, A 7^ 0, 



(2) 



where hvp ^ \/?,ta/2 and a being the C-C bond length. The 
hopping amplitude between the nearest neighboring carbon 
atoms is t ~ 2.8 eV. Due to paramagnetic nature of the non- 
interacting system (i.e. U = 0), the occupation numbers and 
energies of the conduction and valence bands are independent 
of spin orientation. We consider the undoped system with pre- 
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FIG. 1. (Color online) Behavior of free particle hole propagator 
x'^(q, 1^). The imaginary part has been magnified to clearly mark 
the onset of particle-hole continuum of free particle-hole excitations. 
In this figure, vertical scales are drawn for qa = 0.7, and vp and h 
are assumed to be unit. 
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FIG. 2. (Color online) Behavior of free particle hole propaga- 
tor x"(ci, i^) when the gap vanishes, i.e. A = 0. The imaginary 
part (blue line) marks the onset of particle-hole continuum of free 
particle-hole excitations. In this figure, vertical scales are drawn for 
qa = 0.7, and vf and h are assumed to be unit. 



cisely one pz electron per carbon atom, and assume the tem- 
perature to be zero. Then we extend the collective mode anal- 
ysis of Ref. mil to the case with A 7^ 0, whereby the single- 
particle spectrum changes from ±?iui?|k| to the one given by 
Eq. dU. This change in the single-particle spectrum, changes 
the borders of the particle-hole continuum. However, as will 
be show here, at the border corresponding to massless spec- 
trum, interesting collective quanta can be made in presence of 
short range Coulomb interactions. 

Using equation of motion it can be shown that a simple 
RPA-like expression governs the eigen-value equation for the 
collective excitations in the triplet and singlet channels of the 
two-band arising from Dirac fermions in presence of the short 
range Coulomb interactions are given by. 



l±;7x°(q,c^) =0, 



(3) 



where the — (+) sign corresponds to triplet (singlet) chan- 



nel III2II . and the sign difference can be traced back to fermion 
anti-commutation relation. Note that although in the above ex- 
pression, x° is the particle-hole (polarization) bubble, the (ap- 
proximate) bosonic operators which satisfy the above equa- 
tion are of a peculiar form, which does not necessarily coin- 
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cide with the precise particle-hole fluctuation form lull. The 
non-interacting susceptibility employed in the above equa- 
tion, is defined by. 



(q,w) = — 2^— — ^-p — • (4) 



For the dispersion relation (|2|i, the susceptibility has been 
calculated by many authors, the analytic form of which 



reads lUlLM]: 



7rft(w|,|q|2 
z;|,|q|2-cj2-4A2/;i2 



arctan 



2A 



2A/n 



(5) 



^ij|,|q|2 - ^2 

To understand the structure of this function, in Fig. [T] we 
have plotted the real and imaginary part of this function versus 
uj for some representative value of q. As can be seen, the 
imaginary part is non-zero for lo > a;A(q), where 



'^A(q) 



w||q|2 



4A2/^^ 



(6) 



corresponds to the lower boundary of PHC of massive Dirac 
fermions. The singularity of the above inter-band response at 
w = WA (q) is of a weak logarithmic form which could give 
rise to a solution for Eq. ^ only in the triplet channel, and 
for extremely large values of U which is unphysical. The 
however possesses another much more interesting simple pole 
structure at a lower energy scale, uj ^ luq (q) given by. 



wo(q) = WF|q| 



(7) 



Note that this equation defines a line in the plane of a; and |q|, 
which corresponds to the lower boundary of the PHC of mass- 
less (A = 0) Dirac fermions and is determined by the char- 
acteristic velocity, vp the underlying Dirac fermions. This 
could be considered as a property of a "parent" massless Dirac 
system which manifests itself as a simple-pole structure when 
a mass parameter, A is turned on in the non-interacting part 
of the Hamiltonian. At the gapless point, where A = 0, the 
two energy scales at wo (q) and cja (q) merge and the resulting 
singularity of x° will be of inverse-square-root form ifioll . 

As can be seen, to search for solutions of Eq. (O, one has 
to look for the intersection of the real part (green line) with 
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the constant horizontal line Tjj^ where the upper (lower) sign 
corresponds to singlet (triplet) channel. As can be very clearly 
seen in the figure, near wo(q) ~ w_F|q|, where the imaginary 
part of x° is zero, there can be poles both in singlet and triplet 
channel. The simple pole structure of the particle-hole prop- 
agator near this energy scale implies that the poles exist for 
every value of the Hubbard U, no matter how small or how 
large it is. To see this more clearly, note that the real part of 
the susceptibility near ojq (q) behaves as. 



Rex°(q,c^) 



|q|A 



TTVpfi^ — wo(q) 



(8) 



The solutions of the eigenvalue equations (O for arbitrary 
value of U define two collective branches in singlet and triplet 
channel as. 



singlet/triplet (^) = ^^f? I 1 T 



nh?Vp 



(9) 



where upper sign corresponds to the singlet channel, and the 
lower sign corresponds to the triplet channel. These solutions 
are available even for arbitrarily small of U. This formula can 
be interpreted as the spin charge separation in the sense that 
the singlet and the spin-1 modes move at different velocities 
given by. 



Vc/s 



VF [It 



UAcP_ 



(10) 



The mere existence and simple pole structure of the at the 
energy scale ljq is a unique consequence of the gap opening 
in the single-particle spectrum of excitations. The presence of 
gap in the single particle spectrum, pushes the continuum of 
free particle-hole pairs to higher energies, and hence the above 
collective excitations which appear around wo(q) = vp\c^ 
line are well separated from the boundary of PHC and there- 
fore will be protected from Landau damping to the incoherent 
background of free particle-hole excitations. 



DISCUSSION 

The opening of a single particle gap in the spectrum of ex- 
citations in graphene, although pushes the lower edge of the 
particle-hole continuum from the energy scale wq to energy 
scale OJA, but still leaves behind a signature of underlying 
massless Dirac fermions in the form of a singular behavior 
for x° on the line defined by wo(q) = v^lql correspond- 
ing to the PHC boundary of the underlying massless Dirac 
fermions. The simple pole left on this line has two conse- 
quences: (i) Due to the sign change of the susceptibility across 
Wo- there will be collective mode solutions at both singlet 
and triplet channels, for arbitrary values of the Hubbard U . 
The explicit form of the solutions are given by Eq. (|9]l. (ii) 
The different velocities for the two modes is reminiscent of 
the situation one encounters in one-dimension, where a faster 



"spin mode" overtakes the slower "charge mode". Albeit the 
difference is that in ID the divergences in particle-hole bub- 
bles are of the characteristic inverse square root type. In the 
case of 2D massive Dirac fermions, where A is finite, the di- 
vergence in the particle-hole propagator is of a simple pole 
form. When the limit A ^ is taken, the simple pole merges 
with the logarithmic singularity at the boundary of PHC, uja, 
and gives rise again to a inverse-square-root behavior for gap- 
less Dirac fermions in 2D |j To| . It is therefore the non-zero 
value of A which provides a solution in the singlet chan- 
nel in addition to the known solution of the triplet channel 
at A = point ifioi [Till . In this sense, it appears that the 
non-zero value of A enhances the separation of the triplet and 
the singlet modes. Therefore the smallest value of the gap 
parameter. A, totally changes the nature of collective excita- 
tions in graphene. This observation may have implications 
for novel approaches to the bosonization of the 2h-1 dimen- 
sional massive Dirac fermions. Note that if the interaction 
employed was of a long-range Coulomb form, the solution in 
the spin-1 channel would be of a linearly dispersing gapped 
form, and the solution in the singlet channel would not exist 
at all. Therefore the two modes discussed here are peculiar 
feature of short range interactions, which due to remarkable 
value of the Hubbard U in graphene, may have relevance to 
physically fabricated samples of gapped graphene. 

Now let us discuss the nature of the phase transition marked 
by vanishing of the velocity of the singlet mode. As can be 
seen in Eq. ( fTOl i, when the Hubbard U is large-enough to sat- 
isfy. 



?7,A 
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(11) 



the velocity of singlet excitations becomes zero, and the low- 
energy sector is exhausted by only triplet excitations of the 
form 



Wtriplot(q) = 2vFq- 



(12) 



Thinking in the spirit of slave-boson approach lllSll . an elec- 
tron (hole) can be assumed to be composed of its spin part, the 
spinon, and the charge part the doublon (holon). Therefore a 
composite object constructed from an electron and a hole can 
be thought of a separate spin-1 composite of two spinons, and 
a charge composite of a doublon and a holon giving rise to 
a spin zero boson. Hence the triplet mode can be interpreted 
as a triplet bound state of two spinons ifloifTlll . and the new 
singlet mode the emergence of which is being facilitated by 
a non-zero gap in the single-particle spectrum, can be inter- 
preted as a composite boson constructed from a doublon and 
a holon the total charge of which is zero, as it should be. In 
this sense, the vanishing of the velocity of singlet mode can 
be associated with infinite enhancement of the effective mass 
of charge bosons. This gives rise to localization of charge 
carriers and hence is a Mott transition [ flil] . Indeed the de- 
creasing trend Uc oc A~^ as a function of the gap parameter 
A in Eq. ( fTTT ) is in qualitative agreement with our previous 



4 



DMFT result for the phase boundary of the Mott insulating 
phase in the ionic Hubbard model ||9|]. Note that this agree- 
ment holds for very small to moderate values of A. For large 
values of A, the phase boundary in DMFT will be given by 
[/ oc A 10]. However, as long as experimental realization 
of gap parameter in graphene is concerned, A will be on the 
scale of few tens of meV, which is quite small in the scale 
of the hopping term t ~ 2.8 eV. Moreover, the fact that the 
only low-energy excitations left in the system when U goes 
beyond Uc is a another evidence that at Uc the systems be- 
comes a Mott insulator Because in the Mott insulator, the 
only possible low-energy modes are spin excitations. 

When U is further increased, it can be seen in Fig.[T]that the 
real part of x° and (triplet channel) will intersect in an- 
other point which is immediately below the edge of the PHC, 
WA(q), and therefore a second triplet mode at slightly higher 
energy than the first one is expected to form, which is antic- 
ipated due to weak logarithmic singularity at wa- However, 
this mode is not likely to be realized as it requires quite a 
large value of U which falls already in the Mott insulating 
phase, where the ground state maybe totally deformed with 
respect to the initial starting point of massive Dirac fermions 
employed in this work. 

It is interesting to note that, although the velocity of the 
underlying Dirac fermions is vp, the velocity of the spin-1 
mode is always more than vp. Specially at the transition to 
Mott insulating phase, the velocity of the spin- 1 mode will be 
2vf, i.e. the triplet mode (if we call it triplon) moves twice 
faster than the underlying Dirac electrons. This prediction can 
be tested not only in gapped graphene system, but also on a 
platform based on cold atoms engineered to mimic a massive 
Dirac theory at nano-Kelvin energy scales. The existence of a 
mode whose velocity can be twice the velocity of underlying 
Dirac fermions may have interesting implications in teleporta- 
tion of quantum information, as well as in the spin-only forms 
of transport. 

The above two modes discussed here exist for small val- 
ues of U too. Therefore it should be possible to check for 
their experimental consequences. The well isolation of the 
energy scale wo from the PHC edge wa makes the gapped 
graphene more interesting for performing neutron scattering 
experiments. In this work we find that the triplet excitation 
exists below the Mott transition. Moreover, the Mott insulat- 
ing phase has its own magnetic excitations. Therefore in the 
gapped graphene we anticipate a triplet excitation over a large 
range of values of the Hubbard parameter U . Such low-energy 



spin-1 branch of excitations will have a characteristic T/v1 
contribution in the specific heat at constant volume, where Vs 
is the velocity of spin modes. Below the Mott transition, there 
will be another T/u^ contribution coming from singlet charge 
modes. The influence of these modes on various properties of 
gapped graphene remains to be investigated. 
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